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Abstract 

I conjecture that only those states of light whose Wigner function is pos- 
itive are real states, and give arguments suggesting that this is not a serious 
restriction. Hence it follows that the Wigner formalism in quantum optics 
is capable of interpretation as a classical wave field with the addition of a 
zeropoint contribution. Thus entanglement between pairs of photons with 
a common origin occurs because the two light signals have amplitudes and 
phases, both below and above the zeropoint intensity level, which are corre- 
lated with each other. 



I. What states of light are real? 

I shall begin with the following 

Conjecture 1 Real states of light have a positive Wigner function. 

Probably the reader does not believe that this conjecture may be true. 
Therefore I shall rewrite it in a more acceptable form: 

Conjecture 2 In all experiments where the conventional quantum interpre- 
tation involves states of light with negative Wigner functions, the interpreta- 
tion would require only positive Wigner functions if all sources of noise, and 
other nonidealities, are taken into account. 

A. The quantum states of light and their Wigner rep- 
resentation 

The "states" which form the basis for the quantum optical description 
of the light field are the Fock states, which are generated by applying the 
creation operators a\ x to the hilbert-space vector 

io>=nio k ,A>, a) 

which represents the vacuum. The whole Fock space is then spanned by the 
set of vectors 

i {n M » = n m = n ^=(<A) nk ' A io k ,A> , (2) 

k,A k,A y^k,A' 

which represents a state having photons of wave number k and polar- 
ization A. The latter index takes either the value 1 or 2. The most general 
pure state is a superposition of these, that is 

I$) = E^(W,a})IW,a}), £M 2 = i. (3) 

The full set of quantum states is obtained by extending the set |<3>) to mixtures 
of the form 

p = Y,\®)p*(®\, o<p*<i, E p * = 1 - 

2 



The Wigner function of this state is defined as 



W p {{a KX \) = Tr [pW({a kiX })\ , (4) 

where {c^a} are a se t °f complex variables, representing the amplitudes of 
the radiation modes, and 

1 r 

W({a KX }) = II -J / exp [C k ,A(4,A - a K\) - C,A«A - « k ,A)J ^^k,A ■ 
k,A ^ J 

(The integration should be performed with respect to the real and imaginary 
parts of every complex variable.) For instance, the Wigner function of the 
vacuum is 

W ({a k ,A}) = II( 2 A0 ex P (-2|«k,A| 2 ) • (5) 

k,A 

In the hilbert-space representation the electric and magnetic fields of the 
radiation are usually expanded in normal modes, the coefficients of the ex- 
pansion being the creation, A , and annihilation, a kA , operators of photons. 
In the Wigner representation these operators become the amplitudes of the 
modes, which are random variables with a distribution given by the Wigner 
function (assumed positive). For instance the expansion of the electric field 
in free space is 



E(r, t) = \ -JT Re ["k,Ae k ,A exp (ik.r - iut)] , u = c|k|, (6) 

k,A ' L 

where H is Planck's constant, ek,A the polarization vector and L 3 the normal- 
ization volume. 

The Wigner function in nonrelativistic quantum mechanicsB plays the 
role of a pseudoprobability distribution; its marginals with respect to posi- 
tion and momentum separately give the quantum probabilities for each of 
these variables, but the function itself is not positive definite. There are 
great difficulties in interpreting the Wigner function as a true probability 
distribution in quantum mechanics. Nevertheless I propose that for the light 
field the Wigner function may be assumed positive for all physically realizable 
states. This assumption will not put too big difficulties for the interpretation 
of actual experiments, as is argued in the next subsection. 
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B. When is the Wigner function positive ? 

The Fock states (0) form a basis appropriate for the solution of the 
Maxwell equations fulfilled by the electromagnetic field. They play the same 
role as, for instance, the functions {sin(nx)} in the solution of the diffusion 
equation in one dimension, 

df = &l 
dt dx 2 ' 

f(x,t) being the density of the diffusing matter. With boundary conditions 
f(0,t) = f(it,t) = 0, the solution is 

f(x, t) = Y^ c n sin(nx) exp(-nH), (7) 

where the coefficients c n may be obtained from the initial condition f(x,0). 
We need the whole set of functions {sin(nx)} with integer n in order to be 
able to express the solution of the problem as a Fourier series (this require- 
ment is called completeness of the basis) , but it is obvious that not all series 
of the form fl7|) may represent real physical states. In particular the func- 
tions {sin(nx)} themselves are not positive definite (except for n — 1) and 
cannot represent densities. Similarly we may assume that the "Fock states" 
do not correspond to physically realizable states (except (J^)) although all of 
them are necessary in order to represent every possible physical state of the 
radiation field. 

Our assumption contrasts with the frequent interpretation of experiments 
in terms of one-photon states. The one-photon state |1) has the Wigner 
function 

W x {a) = (4|a| 2 - l)Wo(a) , 

where Wo (a) is the vacuum Wigner function for a single mode (|5|) . Taking 
into account the complete set of modes, including the "unoccupied" ones, 
the Wigner function may be written as 

Wi({a k , A }) = (4|« kiA | 2 - 1) W ({a k ,A}). 

Of course this is not always positive. 

Actually the experimental situations!!! in which something like a one- 
photon state has been reported, must necessarily involve the observations of 
wave packets rather than single-mode signals. Indeed, the latter, which fill 
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the whole of space and time, are not at all physical objects. A wave packet 
has the hilbert-space representation 

K(x)) = Ec k ,A^ k - x <Alo), (8) 

where {Cica} are a se ^ °f random, but not independent variables satisfying 

EIC M l 2 = i, 

k,A 

and which, furthermore, are nonzero only for a set of vectors k falling within 
a small ellipsoidal region centred at k'. 

The Wigner function of (Q) is not positive. However, it is necessary to 
bear in mind that there is no way of controlling the moment at which such a 
packet is emitted, in the atomic-cascade situation used for such experiments, 
and this may be taken into account by forming an appropriate mixture of 
such wave-packet states. We have been able to show that the one-photon 
state becomes a mixture having a positive Wigner function.!'! The proof 
will not be reproduced here, but a hint may be got by realizing that the 
mixed state with density operator 

P = ^|1)(1| + ^|0)(0| 
has a positive Wigner function, namely 

Wi(K',a'}) = 2K, a | 2 w (W,a'}) • 

Incidentally, with our approach it is possible to explain, in terms of pure 
waves, one of the most tramatic instances of " corpuscular behaviour" of light, 
the anticorrelation after a beam splittei@-.@ But the wave-particle duality of 
light will not be treated further in this article. 

In summary I challenge the current, although rarely explicit, assumption 
that all density operators represent physically realizable states. In contrast 
I propose that only those states having a positive Wigner function may be 
realizable. On the other hand one may ask whether there are real states such 
that their (positive) probability distribution, W, is not a Wigner function, 
that is no density operator p exists leading to W via eq.(JD . Such states might 
violate the Heisenberg uncertainty inequalities which, in our approach, derive 
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from the existence of a minimal noise which may be controlled in experiments. 
Thus I think that the answer is in the negative, in the sense that the amount 
of noise in the preparation procedure cannot be reduced below the Heisenberg 
inequalities limit. Indeed, it will be usually well above that limit. In any 
case the question will not be studied here further.! 

C. Quantum pure states which are real 

Pure states, of the form (^) , rarely have a positive Wigner function. 
Indeed it is known that if the Wigner function is positive, it is gaussian. The 
proof was given in! for a single mode and generalized to many modes in.! In 
contrast, no general rule is known in the case of quantum mixed states. 

In the rest of this article we shall treat only the restricted, but important, 
kind of real states having a gaussian Wigner function. For simplicity we 
shall consider the ideal (i. e. unphysical) situation where only a single 
mode of the field contains photons, so that the set {^a} contains a single 
member, and the number states are designated simply |0), |1), |2) . . .. The 
full Wigner function will be given by the product of the single-mode function 
with W (a^x) for each of the "unoccupied" modes. The generalization to 
many modes should not be difficult. 

The most general single-mode gaussian Wigner function may be written 
in terms of two real parameters, A and B, and a complex one, a, as follows 



W(a) = - exp {-A(Rea - Rea) 2 - B (Ima - Ima) 2 } . (9) 

If AB = 4, W(a) corresponds to a quantum pure state and if AB > 4 to a 
mixed state. Values such that AB < 4 do not provide a Wigner function, 
that is no density operator exists whence W(a) may obtained via (|). 
Among pure states the case A = B, that is 

W a (a) = (2/7r)exp(-2|a-a| 2 ) , (10) 

is called coherent state, which is an idealized form of the continuous-wave 
laser if a = (and represents the vacuum (|5]) if a ^ 0). The case A ^ B, 
that is 

W a , s {a) = (2/tt) exp[-2e 2s (Rea-Rea) 2 -2e- 2s (Ima-Ima) 2 ], 
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with s real is called squeezed state. The hilbert-space representations of these 
two states are, respectively, 

| a) = exp (ab) — a*a) |0) = V] — = exp ( — |a 2 |/2J \n), (11) 



and 



t, s) = exp {sia} 2 — a 2 )) \a). 



(12) 



States ( PH ) and ( |T2"D are the only single-mode quantum pure states which 
are real according to our criterion. 



D. Stochastic optics 

If all real states of light have a positive Wigner function, we may inter- 
pret that function as an actual probability distribution of the amplitudes of 
the radiation modes. Thus quantum optics becomes a disguised stochastic 
theory, where the states of light are probability distributions defined on the 
set of possible realizations of the electromagnetic field. We propose the name 
stochastic opticM for the stochastic interpretation of quantum optics derived 
from the Wigner function. From another point of view, the stochastic inter- 
pretation provides an explicit hidden variables theory where the amplitudes 
of the electromagnetic field are the "hidden" variables!. 

The most dramatic consequence of stochastic optics is that the vacuum is 
no longer empty, but filled with a random electromagnetic radiation having 
an energy ^tku per radiation mode, on the average, as is shown in eq.(||). That 
radiation corresponds precisely to the additional term introduced by Max 
Planck in his second radiation law (see e.g.t3 ). The picture that emerges is 
that space contains a random background of electromagnetic waves providing 
what we shall call a zeropoint field (ZPF). 

The question whether that theory should be named classical is a matter 
of taste. It might be said classical in the sense that it is a pure wave theory, 
essentially the same developed during the XIX Century. Photons are just 
wavepackets, usually localized in the form of needles of radiation (see eq.(|8|)), 
superimposed to the ZPF. Nevertheless the theory departs from classical 
optics in the assumption that there exists a fundamental noise, the ZPF, 
which cannot be eliminated even at zero Kelvin. I prefer to remain closer 
to the current nomenclature and say that stochastic optics is not a classical 
theory. 
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II. What states of light are classical? 



A. Zeropoint field and detection theory 

The ZPF is usually not directly observable, although indirectly it may 
produce observable effects, as explained below. Observable signals consist of 
additional radiation on top of the "sea" of ZPF. Crucial for the stochastic 
interpretation is the assumption that the ZPF has precisely the same nature 
as signals. Therefore explaining why the ZPF is not directly observable, e. g. 
by firing photon detectors, is a non-trivial problem which will not be studied 
here (see, e. g.,0 ). In the following we state the problem more precisely. 

The hilbert-space formalism of detection is based on normal ordering, 
that is putting creation operators to the left and annihilation operators to the 
right. Let us begin with the ideal case where the radiation field is represented 
by a single mode with amplitude a. Then the detection probability per unit 
time is 



P q oc Tr 



pa^a 



-Tr 
2 



p ja^a + do) — 1 j 

= Jw{{a)[\a\ 2 -\]d 2 a={\a\ 2 - l -) w =n (13) 

where p is the state of the radiation, W(oi) the corresponding Wigner func- 
tion, and n is the "mean number of photons" in the state. The first equality 
derives from the use of the commutation relations and the second is the pas- 
sage to the Wigner representation. The symbol {) w means the average of the 
quantity inside weighted with the Wigner function of the state. 

In the general case of many modes and working in the Heisenberg picture, 
it is straightforward to get, for a point-like detector, 111 

P«(r, t)<xj W({a k }) [1 (r, t; {a k }) - I ] d 2N a k = ( J - J ) , (14) 

where 

I(r,t; {« k }) = ce |E(r,t; {a k })| 2 (15) 

is the intensity for a realization of the field at the position and time (r,t). 
W({a]f_}) is the Wigner function of the initial state, N is the number of modes 
(we should take the limit N — > oo at some appropriate moment) and Iq is 
the mean intensity of the ZPF. (We use the " caligraphic" I for the (random 
variable) intensity in order to distinguish it from the (nonrandom) average 
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Iq.) Writing the electric field, E, in terms of the initial amplitudes of the 
normal modes, {«k}, is usually straightforward although lengthy. It should 
be made for every particular experiment (almost all performed experiments 
with light produced by parametric down conversion have been studied with 
the Wigner function formalism in the series of articles@-ll!). 

Eq . may be interpreted as stating that the detector has a threshold 
so that it only detects the part of the field which is above the average ZPF, 
that is the detector removes the ZPF. The quantum rule (14) is just to 
subtract the mean, a formal procedure which cannot be physical because it 
gives rise to "negative probabilities". The problem is not the huge value of 
the zeropoint energy (the ZPF intensity is about 10 5 w/cm 2 in the visible 
range), because the threshold intensity Iq cancels precisely that intensity. 
The problem lies in the fluctuations of the intensity. For the weak light 
signals of typical quantum-optical experiments the fluctuations of I may be 
such that X < Iq. This problem may be solved as discussed elsewhere ,0 but 
it will not be considered here further. 



B. Classical and nonclassical states of light 

The fact that nonclassicality derives from the existence of the ZPF pro- 
vides a criterion to classify the states of the radiation field. We shall call clas- 
sical (nonclassical) those states where the ZPF is irrelevant (relevant). More 
specifically we will consider a state as classical if the total field E(r,t) can 
be decomposed into two independent parts, E (r,t) and Ex(r,t) representing 
ZPF and signal respectively .□ When this is the case, all optical phenomena 
are associated with the signal alone, and the ZPF may be ignored altogether, 
as is the situation in classical optics. In particular detectors remove precisely 
the ZPF, see (|l4j). The independence of the fields implies that the corre- 
sponding amplitudes, {a^ A } and {aj. A }, are independent random variables. 
If the probability densities of these are iy ({ a k a}) an d ^({^k a})> then the 
density of the total field will be their convolution, that is 

W({a kiX }) = J W 1 ({(3 KX })W ({a KX - (3 KX })d 2N p. 

It is well known that the Wigner function of a state of light may be 
obtained by means of the convolution of the (Glauber) P-function and the 
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Wigner function of the vacuum state when the former exists, that is 

W({a k>x }) = J P({(3 KX })W ({a KX - f3 KX })d 2N (3. (16) 

If we identify Wi({ak,A}) with P({«k,A}) it is clear that the decomposition 
of the stochastic field E(r,i) into two independent parts requires that a P- 
function exists which is positive definite. So, to summarize, the existence of 
a positive P-function is a necessary and sufficient condition for being able 
to decompose the field into independent signal and ZPF parts. This is our 
criterion for a state of light to be classical and it agrees precisely with the 
standard quantum-optical definition of classical state. 

There is only one classical "pure" state, namely the coherent state (the 
vacuum state being a particular case). In the ideal situation of a single mode 
its P-function is 

P(ot) = S 2 (a — a) , 

S 2 being the two-dimensional Dirac's delta. This leads to the interpretation 
that the coherent state represents a deterministic (non-random) signal su- 
perimposed to the ZPF. The Wigner function of the state is given in eq.(|10]) . 

An important kind of classical "mixed" state is chaotic light, whose P- 
function is 

p(M) = n — ex p P-^) , (17) 

and the Wigner function 

p(M) = II ro 2 + ^ exp f-^zr) • ( 18 ) 

^ 7T (2n kiA + 1) \ 2n KX + 1J 

The quantity n^ x is, in quantum language, the mean number of photons in 
the mode k,A, as may be easily proved from the final equality in ([13]) . A 
particular case of chaotic is thermal light, where the dependence of n on the 
frequency uo = |kc| is given by Planck's law. 

III. Entanglement is correlation between quantum fluc- 
tuations 

We are in a position to justify the essential conclusion of the present 
paper, stated in the title of this section. But I shall not attempt here a general 
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characterization of entangled states in the Wigner function formalism, which 
is a rather formidable problem. It is closely related to that of determining 
what density operators correspond to entangled states, a still open problem 
actively investigated in the last few years by workers in quantum information. 
As said above I shall restrict attention to gaussian Wigner functions. 

A. Gaussian Wigner functions involving two radiation 
modes 

We are interested in studying a real state of light consisting of radiation 
in two separated regions of space (see section 1 for the meaning of real). 
Light in every region will contain many modes but, for the sake of clarity, we 
shall study an example involving only two modes with amplitudes labelled 
and /3 k , y, assuming that in all other modes we have just ZPF. A more 
physical state would correspond to having two wave packets, the first (second) 
containing many modes with wavevectors close to k (k ) . For simplicity we 
will remove the subindices and label the amplitudes a and j3 in the following. 

We shall consider a classical two-modes state whose marginal for every 
mode corresponds to chaotic light (|T7|). Its P-function is gaussian 

ab 2 

PvifaP) = exp {-a |«| 2 - b\(3\ 2 + c(a(3* + (3a*)\ , (19) 

a, b and c being three real numbers fulfilling a,b > 0, c 2 < ab. The latter 
condition guarantees that P goes to zero at infinite, a necessary condition for 
normalizability. In order to simplify the argument I will consider a slightly 
less general state by assuming a = b, although the generalization to a ^ b 
is trivial. The intensity of the radiation, above the ZPF, is measured by the 
integral ( see (p~3|) ) 



2n 



/ P 12 (a,P) (\a\ 2 + \(3\ 2 ) d 2 ad 2 (3 = - T —^, (20) 

J CL C 



where n represents the mean number of photons per mode. 

The Wigner function of this state is given by the convolution with the 
vacuum Wigner function ( see (|16|).) We shall write it in terms of n and a 
correlation parameter, x, defined by 

2c . , 

x= 5 21 

2a + a 2 - c 2 K ' 
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We get 



W 12 {a, f3) = exp {—A [\a\ 2 + \P\ 2 - x (a(3* + (3a*)] },A = 



(2n+l)(l 
(22) 

This Wigner function represents a classical state if |c| < a, which implies a 
restriction on the range of values of x. In fact, from (^) and ( PH ) we obtain 

An 



An 2 -x 2 (2n + l) 2 ' 



The condition a > leads to 

classical : n > 0, \x \ < . (23) 

~ " 1 2n + l v ; 

A transparent interpretation of these results emerges in our approach. 
Remembering eq.flTED , we see that (22) is the probability distribution of the 



sum of two random variables representing the signal and the ZPF, respec- 
tively. In the signal, whose distribution is (0) , there is correlation between 
the two modes. In contrast in the ZPF, whose distribution is (J5|) , all modes 
are uncorrelated (the distribution is a product of single-mode terms). In 
summary, we have a correlation involving the signal but not the ZPF. This 
should be called a classical correlation. 

The constraints for the state (^2|) to be real are weaker than (0), namely 

real : n > 0, \x\ < 1. 

(Another condition for the state to be real is that W may be obtained from 
a density operator via (|j), which is true in this case but will not be proved 
here). On the other hand the marginal of 



Wi (a) = — exp 



2 a 



7T (2n + 1) { (2n + 1) J 

represents a classical state. I shall call entangled any state of two modes 
which is not classical but has classical marginals, that is 

2n 

entangled : n > 0, - — -j— j- < \x\ < 1. (24) 
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(I propose this as a sufficient condition, not excluding the possibility of en- 
tangled states with nonclassical marginals). We see that such states involve a 
correlation between the two modes which is larger than the correlation of any 
classical state (given by (p3|)). This is because not only the signal but also 
the ZPF or "quantum vacuum fluctuations" are correlated in the entangled 
state. 

"Two-photon entanglement" similar to this one (although involving 
many modes) occurs in the process of parametric down conversion . In fact, 
light produced in that way consists of two separated beams whose Wigner 
functions are gaussian but the full state is not classical. That is, every beam 
alone consists of chaotic light of the type flT8|), but the two beams are entan- 
gled§-.§ 

The reader may realize that entanglement becomes less relevant when 
the intensity n (the mean number of photons) is large because the interval 
of x in fl2l|) becomes narrow. The opposite is true if n is small. This is to be 
expected because large n corresponds to the classical limit. 

B. Entanglement and Bell's inequalities 

There is no agreement about the definition of entanglement for mixed 
quantum states. A fashionable criterion is the violation of a Bell inequality. 
In actual experiments the inequalities tested have never been genuine Bell 
inequalities, derived using only general properties of local hidden variables, 
but inequalities involving auxiliary assumptions. The violation of one of such 
inequalities does not imply the refutation of local realising a fact qualified 
as "existence of loopholes" . Nevertheless we may consider the experiments as 
valid tests of entanglement. In practice any test involves measuring a coin- 
cidence detection rate as a function of some controllable angular parameter, 
0. The inequalities are violated if the measured coincidence rate, i?i 2 , is of 
the form 

R12 = const, x (1 + Vcos<p) , (25) 
with the visibility or contrast, V, greater than some limit, usually 0.71. 

The most frequent tests involve polarization correlation. In order to study 
polarization we need to take at and (3 as (complex) two-dimensional vectors 
or, what is equivalent, to double the number of modes. Introducing the 
polarization vectors, e k A , see (|5]) , the vector amplitude a might be written 

a =a x e K1 + a y e K2 , |e*| 2 = \a x \ 2 + \a y \ 2 , 
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and similar for (3. The Wigner function will be, instead of 



W 12 {a,P) = (1 X P A \ x P {-A[\(x\ 2 + \(3\ 2 - x (a ■ (3* + (3 ■ a*)]} , 



71 ' 



A = (2n + 1) (1 - x 2 ) ■ (26) 

We shall consider that the amplitudes a and /3 correspond to two polarized 
light beams arriving at two polarization analizers at angles <p 1 and </> 2 , re- 
spectively. In what follows we may ignore the ZPF in all modes except those 
included in ( p6|) . Consequently we assume that the amplitudes emerging from 
the polarizers are given by Malus law, that is 

A = (a • m) ,fi = {(3- u 2 ) , 

the vector u x having components (cos0 1 , sin^) and similar for u 2 . Note 
that the scalar amplitudes A and n correspond to modes with polarization in 
the directions of U! and u 2 , respectively. As said above we ignore the modes 
with polarization perpendicular to u x or u 2 , which contain just ZPF. 

The coincidence detection rate in two detectors placed after the polarizers 
may be calculated by a straightforward generalization of flT3| ) , namely 

R12 oc J W l2 (a,f3) (jA| 2 - -) - d 2 a x d 2 a y d 2 (3 x d 2 (3 y . 



The integration is trivial using eq.(Ea) and we get 



1 / 1\ 2 

R 12 oc n 2 + - (n + - x 2 [1 + cos (20 x - 20 2 )] 



2 



The visibility is 



V 



'" + |) x<2 



n 2 + \(n + \) 2 x* 

which is greater than | for entangled states, fulfilling condition ( |2~4] ) , but 
it is smaller than | for classical states, where condition (|23| ) holds true. In 
particular, the limit V = 0.71 may be surpassed if n < 0.42 ( but this is 
specific for the gaussian Wigner functions studied in this paper.) 

We see that high visibility is only possible with weak signals. I conjecture 
that the signal weakness, combined with the necessity of removing efficiently 
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the ZPF, gives rise to the difficulties for performing " loophole- free" tests of 
Bell's inequalities. This problem will be studied elsewhere. 
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